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Abstract 

^ In a superstring framework, the free energy density J-" can be determined unambiguously 

(-| at the full string level once supersymmetry is spontaneously broken via geometrical fluxes. 

We show explicitly that only the moduli associated to the supersymmetry breaking may 

give relevant contributions. All other spectator moduli /// give exponentially suppressed 

^ contributions for relatively small (as compared to the string scale) temperature T and su- 

00 persymmetry breaking scale M. More concisely, for jii > T and M, J-" takes the form 
00 

^ J-(T, M; = ^(T, M) + o[ exp (-^) , exp [-%)'■ 

O 

^ We study the cosmological regime where T and M are below the Hagedorn temperature scale 
Th- In this regime, J-" remains finite for any values of the spectator moduli yU/. We investigate 
. ^ extensively the case of one spectator modulus Hd corresponding to Rd, the radius-modulus 
^ field of an internal compactified dimension. We show that its thermal effective potential 
V{T, M; jj) — T{T, M; jj) admits five phases, each of which can be described by a distinct 
but different effective field theory. For late cosmological times, the Universe is attracted to a 
"Radiation-like evolution" with M{t) oc T(t) oc l/a{t) oc t"^/*^. The spectator modulus ^(t) 
is stabilized either to the stringy enhanced symmetry point where i^^ = 1, or fixed at an 
arbitrary constant //q > T,M. For arbitrary boundary conditions at some initial time, tE, 
li{t) may pass through more than one effective field theory phase before its final attraction. 
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1 Introduction 



String theory provides a framework to obtain a sensible theoretical description of the cos- 
mological evolution of our Universe. Nowadays, it is the only known framework in which the 
quantum gravity effects are under control [l], at least for certain physically relevant cases. 
Following the stringy cosmological approach developed recently in Refs (2]-^, the classical 
string vacuum is taken to be supersymmetric with a fixed amount of supersymmetries defined 
in flat space-time. 

This initial choice does not give rise to any cosmological evolution. In the presence of 
supersymmetry, the quantum corrections to the gravitational background would lead to a flat 
space-time, or would modify it at most to Anti-de Sitter, domain walls or gravitational wave 
backgrounds respecting a time-like or light-like killing symmetry. The above cosmological 
obstructions are however physically irrelevant for two fundamental reasons: 

• Firstly, supersymmetry is broken in the real world, (at least spontaneously and not 
explicitly), at a characteristic supersymmetry breaking scale M . 

• Secondly, in the case of thermal cosmologies, the supersymmetry is effectively (spon- 
taneously) broken at the temperature scale T. 

Both the M and T supersymmetry breaking scales induce at the quantum level a non- 
trivial free energy density J-'{T,M), which plays the role of an effective thermal poten- 
tial V{T,M) = J-'{T,M) that modifies the gravitational and field equations, giving rise 
to non-trivial cosmological solutions, as has been explicitly shown in Refs [2}[3|[5]. Both 
the supersymmetry breaking and finite temperature phenomena can be implemented in the 
framework of superstrings [6]-[8] by introducing non-trivial "fluxes" in the initially supersym- 
metric vacua. Furthermore, in the case where supersymmetry is spontaneously broken by 
"geometrical fluxes" [9 ,10 , the free energy J-'{T, M) is under control and is calculable at the 
full string level, free of any infrared and ultraviolet ambiguities (2|[3]. This is true, provided 
T and M are below a critical value close to the string mass scale, the so-called Hagedorn 



temperature Th [6j]8 11 . In the framework of stringy-thermal cosmologies, T ~ Th corre- 
sponds to very early times when we are facing non-trivial stringy singularities indicating a 
non-trivial phase transition at high temperatures |6|[8| [l2|[l3|[l5] at time tn- In the literature, 
there are many speculative proposals concerning the nature of this transition [6|[8| [l2}jl6] . 
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A way to bypass the Hagedorn transition ambiguities was proposed in Ref. [2]. It con- 
sists of assuming the emergence of {d — 1) large space-hke directions for times t ^ tff, 
describing the {d — l)-dimensional space of the Universe, and possibly some internal space 
directions of an intermediate size characterizing the scale M of the spontaneous breaking 
of supersymmetry via geometrical fluxes |9 , 10 . Within these assumptions, the ambiguities 



of the "Hagedorn transition exit at ts" can be parameterized, for t > t^, in terms of ini- 
tial boundary condition data at t^; > tn- In this way, the intermediate cosmological era 
tE < t\Y, i.e. after the "Hagedorn transition exit" and before the electroweak symmetry 
breaking phase transition at t^, was extensively studied in Ref. j4] in the case of c? = 4. An 
output of the present analysis is that the cosmological "radiation-like" evolution found in 



Refs [2|[3|[5 17 generalizes to a "Radiation-like Dominated Solution" (RDS) in c?- dimensional 
space-time, 

RDS^ : M(t) oc (t) oc l/a{t) oc t"^/^ for ts < t < tw, (1-1) 

and is unique at late times in certain physically relevant supersymmetry breaking schemes. 
As a necessary and sufficient consistency requirement, we note that in this intermediate 
cosmological regime Th ^ T, M ^ Tw, the smallness of the space-time curvature scales, 
= (a/a)^ and H, the dilaton and the evolving radii scales, ^^n, (t>dii,{Ri/Ri)^, Ri/Ri, are 
guaranteed to be small (< 0{T'^, M'^)), thanks to the "attractor mechanism" towards the 
RDS in late cosmological times. In particular, they are all decreasing at late cosmological 
times and so our quasi-static approximation becomes better and better as time passes. In 
addition, our perturbative approximation becomes better and better as time progress due 
to the falling of the dilaton. 

We point out that the evolution is radiation-like in the sense that the external space- 
time's evolution is identical to that of a radiation dominated universe. However, in our 
case the expansion is driven not only by radiation, but also by the coherent motion of the 
supersymmetry breaking modulus M. We would also like to stress that a key result of |4] 
was that the RDS^ is actually an attractor of the dynamics. The important consequence of 
this attractor, is to wash out the dependence of the cosmological evolution on the choice of 
initial boundary conditions, which were used to parameterize our ignorance of the physics 
involved in the Hagedorn transition. In this paper we show that the attractor naturally 
extends to the higher-dimensional RDS'^. 
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Although this analysis was done in the framework of initial vacua with A/4 = 2 super- 
symmetry, the claim is that it will still be valid in more realistic models with initial = 1 
18 . We would like to stress here that the limitation t <tw the infrared 



supersymmetry 



regime follows from the appearance in the low energy effective field theory of a new scale, 
namely the "infrared renormalization group invariant transmutation scale Q" , at which the 
supersymmetric standard model Higgs (mass)^ becomes negative, (no-scale radiative break- 
ing of SU{2) X U{1) f/(l)em |l9l|20)). Q is irrelevant as long as M,T > Q; however 
it becomes relevant and stops the M{t) evolution when T ^ Q at t ~ tp^/ z.e. when the 
electroweak breaking phase transition takes place. Although the physics for t ^ tw is of 
main importance in particle physics and in inflationary cosmology at tw, it will not be ex- 
amined in this work. The main reason for us is its strong dependence on the initial vacuum 
data which screens interesting universality properties. We therefore work in the intermediate 
cosmological era tE < t < tw or Th ^ T ^ Q, i.e. after the Hagedorn phase transition 
and before the electroweak one. In this regime the transmutation scale Q can be consistently 
neglected and, furthermore, the Hagedorn transition ambiguities are taken into account in 
terms of initial boundary conditions (IBC) after the "Hagedorn transition exit" at tE- This 
scenario gives a dynamical explanation of the smallness of the supersymmetry breaking scale 
as compared to the string or Planck scales. Indeed, extrapolating the RDS'^^'^ up to the low 
energy regime where T = 0{1 TeV) one finds, (thanks to the attractor mechanism), that the 
natural value of the supersymmetry breaking scale M{t) is naturally small and around the 
electroweak phase transition, independently of its initial value at early cosmological times. 

The only known supersymmetry breaking mechanism that can be unambiguously adapted 
at the string perturbative level, is the one we consider here where supersymmetry is spon- 
taneously broken via "geometrical fluxes" . This choice implies the existence of at least one 
relatively large compact dimension (the one which is associated to supersymmetry break- 
ing). This is by far not in contradiction with experimental results both in particle physics 
and cosmology. For several years, the possibility of "large extra dimensions" has attracted 
the attention of the particle physics community; the future data analysis at the LHC and 
elsewhere includes searches for signals indicating the existence of large extra dimensions at 
scales 1/R ~ 0{1 TeV), which is the characteristic prediction of supersymmetry breaking 
via geometrical fluxes. 
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Many other choices for supersymmetry breaking exist. However in most cases it is not 
well known, in our days, the precise stringy corrections to all orders in a' that are necessary 
in order to study the intermediate cosmological regime. In all other supersymmetry breaking 
mechanisms we are forced to work in the effective supergravity framework. Our results may 
help to make more general stringy approaches possible in the future. There are indications 
that our results can be converted to other string vacua, utilizing string/ string and M-theory 
dualities. In particular the geometrical fluxes are mapped to other types of fluxes, for 
instance the three form R-R and NS-NS fluxes in type IIB - orientifolds. In this respect, 
the exact stringy approach in the intermediate cosmological regime gives us profound (non- 
perturbative) information via M-theory and string dualities. 

In addition to generalizing the results of |4] to arbitrary dimension, we analyze the 
time behavior of the spectator moduli not participating in the breaking of supersymmetry. 
Following Refs [2|[3|[5], one can show that only the supersymmetry breaking moduli M{t) 
and T{t) can give a relevant contribution to the free energy density J-". Intuitively, all other 
moduli /i/ are either attracted and stabilized to the "stringy" extended gauge symmetry 
points, close to the string scale fij ~ Mgtring, or are effectively frozen to an arbitrary value 
such that fij ^ T and M, giving rise to exponentially suppressed contributions: 



One point of the present paper is to explicitly verify this intuition for the moduli coming 
from the spectator tori. Indeed, the supersymmetry breaking moduli generate a non-trivial 
potential for the spectator moduli and freeze them as expected. Considering the effect on 
a single spectator modulus (/i oc 1/Rd) in the case of the heterotic string, the thermal 
effective potential V{T, M; fi) admits five distinct phases, each of which can be described 
by a different effective field theory. The interesting result is that by using a string theory 
framework and working at the full string level, we are able to link together, within a single 
framework, the different effective field theories. Furthermore, the main result of this paper 
is to derive at the string perturbative level, (however exact in a'), the full string free energy 
J-'(T, M; fij) as a functional of the supersymmetry breaking moduli T, M, the string coupling 
constant modulus and "spectator moduli" fij, for a certain class of string vacua where 
the spontaneous breaking of supersymmetry is induced by geometrical fluxes. 

The form of the potential is sketched in Fig. [T] and the phases are summarized as follows 




(1.2) 
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Fi gure 1: Qualitative shape of the effective potential V versus InRd, with T , M (and the dilaton in Einstein 
frame) fixed. When Ini?^ varies, five distinct phases can arise in the heterotic case, i) The Higgs phase I, 
ii) the flat potential phases II, IV and Hi) the higher- dimensional phases 111, V. The phases IV and V are 
T-dual to II and III. In the type II case, the phase I does not exist and the plateaux II and IV are connected. 
The plot is valid for T > M , for T < M one simply replaces Rq with Rg. 



I. Higgs phase: With 



Rd — 



1 



1 A/ 1 

< -s- and/or — . 

Kq Kg 



This phase contains the stringy extended symmetry point at the self-dual point Rd = 1. 
The appropriate effective field theory description of this phase is in terms of a rf- dimensional 
theory of gravity coupled to an SU{2) gauge theory. R^ is dynamically stabilized at this 
point. Such a notion of moduli stabilization, has been studied in the literature before 
Here we demonstrate such moduli stabilization in the context of the heterotic 



13 21,22 



superstring by an explicit computation of the effective potential. 

1 1 1 

II. Flat potential phase: With — and — < Rd — — < Rq and -Rg. 

Rq Rg Rd 

Here, the appropriate effective field theory description is in terms of a (i-dimensional theory 
of gravity coupled to an f/(l) gauge field. 

III. Higher-dimensional phase: With Rq and/or Rg < Rd. 

For macroscopic values of the spectator modulus Rd, the appropriate effective field theory 
description is the {d + l)-dimensional theory of gravity. The modulus Rd becomes the gdd 
component of the metric, Qdd = {'^'^RdY and the evolution is attracted to that of an RDS in 
d + 1 dimensions. 

One may also consider the case in which the radius Rd is still internal. For large enough 
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values of Rd so that we can neglect terms of order {Ro/RdY and {Rg/RdY, the evolution is 

attracted to an RDS'^"''^ for a long period of time. However, at late times Rq and Rg always 

catch Rd and the solution is ultimately attracted to the RDS'^ of phase II. 

1 1 1 

IV. Dual fiat potential pliase: With — and — < Rd < Rq and Rg. 

Rq Rg Rd 

The effective theory description is T-dual to that of phase II. The light degrees of freedom 
are the winding modes instead of the Kaluza-Klein momenta of phase II. 

V. Dual higher-dimensional phase: With Rq and/or Rg < —. 

Rd 

This phase is T-dual to phase III. Its properties are derived from the ones of phase III under 
the replacement Rd ^ 1/Rd and phase II — )■ IV. In particular, the dual effective field theory 
is a (c? + l)-dimensional theory of gravity, with gdd = (Svr/i?^)^. 

The above different phases of a common string setting cannot be described in the context 
of a single field theory. This is due to the necessary presence of the string winding modes. 
In a field theory framework, only the phases II and III (or IV and V) can be described by a 
common field theory. The winding modes are particularly important for the stabilization of 
the modulus in phase I at the extended symmetry point, and furthermore for the description 
of the T-dual phases IV and V. In contrast to field theory, string theory naturally interpolates 
between these various phases, due to the generation of an effective potential in the presence 
of temperature and spontaneous supersymmetry breaking. 



For each phase, there exists an RDS as in (1.1). We show that these solutions are stable 
against small perturbations and that for arbitrary IBC close to an RDS, the cosmological 
evolution is attracted to this RDS. In |4|, the spectator moduli were taken to be frozen and 
it was shown that under this hypothesis the RDS is a global attractor. Taking into account 
the existence of gravitational friction for an expanding universe, we expect in the present 
work with dynamical spectator moduli such as Rd, the results of j4j to generalize so that the 
evolution is always attracted to the RDS of one of the five phases. 

In the case of type II strings, the phase I does not exist perturbatively, so that the 
phases II and IV combine into a single plateau. This is due to the lack of massless states 
necessary to enhance the U{1) to SU{2). However, by heterotic-type II duality, we expect 
such an enhancement to exist non-perturbatively, so that all five phases should exist at the 
non-perturbative level. The effects correspond to the addition of branes whose separation is 
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governed by the spectator moduli. 

The organization of the paper is as follows. In section 2, we discuss in more details the 
specific setup analyzed in this paper. The thermal effective potential is given and it is shown 
to have the five phases discussed above. In section 3, we show that an RDS solution exists in 
each phase. In addition, we show the stability of the solutions against small perturbations. 
In section 4, we briefly discuss the role of non-perturbative objects in the type II string 
theory. In section 5, we summarize our results and discuss further avenues of research. 
In appendix A, the thermal partition functions for the heterotic and type II strings are 
presented, together with their asymptotic properties to be used in the different phases. The 
gravity and field equations are given for each effective field theory phase in appendix B. 



2 Effective thermal potential in superstrings with 
spontaneously broken supersymmetry 

In the presence of temperature, the one-loop partition function of both the heterotic and 
type II superstrings is non-vanishing and yields the one-loop effective potential at finite tem- 
perature. In addition, spontaneous supersymmetry breaking is induced by the presence of 
geometrical fluxes l2l[3] along the internal cycles of the background manifold. We introduce 



these fluxes via the generalization to the context of string theory of Scherk-Schwarz 23 



compactifications in field theory 24 , 25 . They induce further contributions to the one-loop 
partition function, which persist even at zero temperature |2||3]. Due to various ways of in- 
troducing the fluxes, there are multiple supersymmetry breaking configurations for the same 
initially supersymmetric background. In [2||3], such one-loop thermal effective potentials 
were derived in the limit of small temperature T and small supersymmetry breaking scale 
M for the heterotic and type II superstrings compactified on and x — orbifolds. 



Z2 

The partition functions were calculated for small but otherwise arbitrary temperature and 
supersymmetry breaking scale, while the remaining moduli were taken to be frozen close to 
the string scale. 

We want to relax the latter hypothesis and examine the behavior of the spectator moduli 
in the presence of temperature and supersymmetry breaking. In appendix A, we compute 



the partition functions for the heterotic (A. 16) and type II (A. 21) cases. The background 
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manifold is of the form S]^ x x (or S]^ x x — x in the orbifold models) , where S]^ 
is the compact Euclidean time circle and the T" torus involves the geometrical fluxes which 
generate the breaking of supersymmetry. The spectator moduli are not participating in 
the breaking of supersymmetry. We take both the temperature and supersymmetry breaking 
scales to be small, while allowing the spectator moduli to remain arbitrary. This enables 
us to study the resulting effects of the effective thermal potential on the moduli. 

For simplicity, we specialize to the following 10-dimensional Euclidean background that 
contains: 

• The Euclidean time direction, with radius Rq which determines the temperature T. 

• The 1, ...,d — 1 directions, which are taken to be very large and form, together with 
the time, a c?-dimensional space-time. 

• The circle S^{Rii), with arbitrary radius. For small R^, S^{Rii) is considered as part of 
the internal compactified space. For macroscopic Rd, however, S^{Rd) becomes part of 
a space-time of dimension d + 1. By macroscopic we mean that we can probe it with 
current experiments. Rd is the only "spectator" radius whose dynamics is taken into 
account. 

• The n = 1 circle involved in the spontaneous breaking of supersymmetry. We take it 
to be along the compact direction 9, with radius Rg. 

• The remaining compact directions, with radii Rd+i, -Rs- They are taken to be fixed 
close to the string scale. (In the orbifold models, the — factor spans the directions 
5, . . . , 8. Its dynamics are consider in the companion paper [18].) 

Utilizing the general expressions associated with the heterotic and type II partition functions 
given in appendix A, we can easily obtain the ones associated to the background chosen 
above, namely 

S\Ro) X T'^-i X S\Rd) xMx S\Rg), (2.1) 

where Ai = T^^'^ or T^"*^ x — . The heterotic (type II) models admit a supersymmetry 
characterized by 16 or 8 (32 or 16) supercharges, which are spontaneously broken by the 
"stringy Scherk-Schwarz compactifications" in the directions 9 and (2|[3]. The scales of 
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supersymmetry breaking M and temperature T are characterized by I/-R9 and I/Rq, respec- 
tively. We take Rq and -Rg to be large, but still much smaller than the radii of the T'^~^ 
torus so that we have the following inequality 

-Ri, -Rd-i > -R05 -Rg > 1- (2-2) 

As long as Rd is smaller than the size of the external space T*^"^, it is more convenient 
to express the effective field theory action S in terms of fields, which have a natural inter- 
pretation in d dimensions. We are interested in isotropic and homogeneous backgrounds. 
More specifically, we take the gauge fields to be pure gauge and the remaining scalar fields 
and the space-time metric to depend only on time. This will have the advantage that af- 
ter such a reduction, the different effective field theories will be describable within a single 
framework. The backgrounds we will consider are non-trivial for the c?- dimensional metric 
g^,^, the (i- dimensional dilaton ^^ii, and the moduli fields. However, since we allow Rd to 
vary arbitrarily in size, it may become of the order of the external T*^"^ radii, so that S^{Rd) 
should be considered as a part oi a {d+ l)-dimensional space-time. In this case the effective 
action S is naturally expressed in terms of redefined fields and space-time metric in d + 1 
dimensions. 

In appendix [B| the dimensional reduction from 10 dimensions to d dimensions is carried 



out explicitly and the resulting action in Einstein frame is given in (B.7). The case we are 
considering here has n = 1, A = 9 — A~d (where A = in the toroidal models and A = 4 
in the orbifold ones) and T> = d. The resulting action is 

s= j d'xv^g (I - lid^r - lid<i>±r - lidcr + , (2.3) 

where we have defined the normalized fields, 



along with 

C:=\nRd, r]:=lnRg. (2.5) 
The remaining moduli are taken to be fixed close to the string scale. 



The source P is a pressure equal to —J-", the free energy density (see Eq. (B.12)). It is 



the opposite of the one-loop effective potential at finite temperature and is related to the 
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one-loop partition function as, 

p ^ gd32</'dil ^ (^2.6) 

Vo,...,d-i 

where Vo,...,(i_i is the d-dimensional Euchdean volume (in string frame) and Z is the partition 
function computed in appendix A. In the next sub-section we shall give the exact form of P 
in terms of a convenient set of variables. 



2.1 Specific form of the effective thermal potential 



For the heterotic case, the partition function is given in Eq. (A. 16) which can be re-written 
in the following convenient form 



Z — j Ri 1 Rg ^ ^ (^generic ^enhanced) 



(2.7) 



90+S9=l 



For the type II case, Eq. (A. 21), there is no contribution 



■'enhanced * 



In the heterotic case, 2^onhancod is generically suppressed except when Rd is close to its 
self-dual point, where an enhancement of the gauge group U{1) — )• SU{2) occurs. The 
contribution ^generic for generic Rd can be written in two equivalent form^ related to one 
another by a Poisson resummation of the momentum lattice index nid of the circle S^{Rd)- 

In the Einstein frame, the temperature T and supersymmetry breaking scale M are 



dressed by the dilaton field (p^n and are given by Eqs (B.IO) and (B.5): 



T 



27rR 



M 







e d---^ 
2nRo, 



2-K 



(2. 



M Ro 



the (piii dependence drops out. The expression for P 



Observe that in the ratio ^ 

T Rg 

gets simplified drastically once it is written in terms of the complex structure modulus z, 



z :=\n 



M 



In 



Rn 



(2.9) 



T Rg 

In terms of the independent variables {T, z, t], (}, the pressure P takes the factorized form 



P{T,z,r]X)^T'p{^,vX), 



(2.10) 



^In the notations of appendix A, these two forms correspond to A — 9 — A — d and A = 8 — A — d, where 
A — ioi the toroidal models and yl = 4 for the orbifold ones. 
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with ( and 1] defined in (2.5). Furthermore, p can be written in terms of functions with 
natural interpretations either in or + 1 dimensions. (In Eq. (A. 16), the first 



case 



corresponds toA = 9 — A — d and the second to A = 8 — A — d.) In the heterotic case, the 
two equivalent forms for p are: 



id), 



+ nv 



.(d), 



+ nT9T\z,V, ICI) +^y^r(^'^' Kl) 



(d), 



riT 



f!f^'\z^r,-\C\)+nyf';f^'\z,r,-\C\) 



(2.11) 



+ nTgT\z,r], \C\) +nvgv\z,r], |C|), 

with the functions defined below. For the type II case, one simply takes ut = nv = 0. Note 
that p is an even function of (, as follows from T-duality Rd — )■ 1/Rd- In this expression, 
riT is the number of massless boson/fermion pairs of states in the originally supersymmetric 
background, for generic Rd- ut is the number of additional ones at the enhanced gauge 
symmetry point. The value of ny is given by the sum over the ut pairs, with each pair 
weighted by a sign. The distribution of signs depends on the specific supersymmetry breaking 
configuration and can yield a negative ny. For the heterotic models we consider, one has 



2^ 



-1<^<1 



flT _ 2 



nv = riT ■ 



(2.12) 



The definitions of the various functions appearing in Eq. (2.11) are given by 



TT 2 



[e2-(2^o + 1)2 + i2kgy 



d+1 5 
2 



4'^(^,r/-ici) = E' 



"Id 



ko,kg 



2Kd±i (27i\md\e''-\^\ ^Je^' 



e^'{2ko + 1)2 + (2A;9)2 



62^(2^0 + 1)2 + (2A;9)2 



#(.,ry,|C|) = (e2KI-l) 



2|CI 



d+l 



e 2 



--{ri-KDdz 



ko,kg 



2K^ 27r(e2|fl - l)e''-\'^\ e^'{2ko + 1)2 + (2^9)2 



e2^(2A;o + 1)2 + (2A;9)2 



d+l 
4 



d 

7T2~ 



ko,k<),rha 



[e2-(2fco + 1)2 + (2^9)2 + 6- 



(2.13) 
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where Ka{z) are modified Bessel functions of the second kind. The remaining functions with 
lower index V are related to those with lower index T by M ^ T duality transformations 
{z ^ -z): 

k'v\z,v- ICI) = e^'-'^^ k!^\-z,v- ICI + ^) 

(2.14) 

flf^'\z,V- ICI) = e''^ f^^^'\-z,^- Id + z) . 

We will first focus on the dynamics of the modulus Rd, and so we consider the behavior 
of —P = —T'^p{z, rj, () at fixed T, z and rj. Note that when the functions g!^^ and gy^ can 
be neglected, the pressure only depends on two quantities, z and r] — \(\. 



2.2 The five heterotic effective field theory phases 

Considering the effect on a single spectator modulus Rd in the case of the heterotic string, 
the thermal effective potential —P admits five distinct phases corresponding to different 
effective field theories. The form of the potential is sketched in Fig. [T] and the phases are 
summarized as follows: 



• I: Higgs phase 



Rd 



1 

Rd 



< -s- and/or — . 

Uq Kg 



(2.15) 



This phase contains the stringy extended symmetry point at the self-dual point -R^ = 1- 
The appropriate effective field theory description is in terms of a (i-dimensional theory 
of gravity coupled to an SU{2) gauge field. We will see that the modulus Rd can be 
stabilized at the self-dual point which turns out to be the minimum of the effective 



thermal potential. Indeed, considering the expression of p in (2.11), the functions gj. 



and / or gy^ are of the same order as f^' and fy' , while fc^'' and ky' are exponentially 
small due to the behavior of the modified Bessel functions Kd+i. In particular, one 

2 

has at the origin ^ = In i^^^ = the following behavior. 



(d) 



p{z,V,C = 0) = {nT + nT)fT\z) + {ny + ny) fy' (z) ■.= p{z). 



id), 



(2.16) 
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This is precisely the form obtained in p], when the dynamics of Rd was ignored i.e. 



Rd was taken to be stabihzed close to the string scale. In Eq. (2.16), the contribution 
of ut and ut is of the same form since both contributions come from massless states 
when ( is at the enhanced gauge symmetry point C = 0- Due to the fact that ht and 
hy are positive, the extremum of —P at C = is always a minimum. 

II: The fiat potential phase 

1 1 1 

— and — < Rd--s-, 

Kq Kg Kd 

and Rd < Ro and Rg. (2.17) 

For this range of the modulus, there exists a description in terms of a ci-dimensional the- 
ory of gravity coupled to a U{1) gauge field. Note that the range of this region grows as 
Rq and Rg increase. Modulo exponentially suppressed terms O [ exp(— ^), exp(— -^) ] , 
the potential for the modulus is flat. We will see that for certain IBC, the modulus Rd 
may be frozen to an arbitrary value on this plateau, due to the gravitational friction 
of the expanding universe. The exponentially suppressed terms are irrelevant in this 
phase and cannot modify this behavior. 



In this range (2.17), the contributions of k!^'' and ky\ as well as g^'' and gy'' are 



exponentially small compared to f!^^ and fy \ so that —P is independent of ( = \nRd. 
The pressure reproduces the result in d dimensions for ht massless boson/fermion 
pairs in the originally supersymmetric model (as opposed to phase I which has ut + nr 
massless pairs). Physically, we are away from the enhanced symmetry point and so 
the previous SU{2) states are no longer massless. More concretely, along the plateau 
we have 

p{z, r^, ^ riT f^\z) + nv flf\z) := p{z) . (2.18) 
Either sign of this quantity is allowed when ny < 0. Indeed, considering the large \z\ 



behavior, (2.18) implies 

-P ~ -T'^nre-' ^ -oo 



-oo 



~ -T'^ nv e'^^ S'^,^ sign(-ny)oo 

->+oo z— >+oo 



(2.19) 



where (and S"^*^ for later use) is a constant. 



^E^. 5"-^E'i. (2.20) 
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and we see that —P may take any value. 



• III: Higher-dimensional phase 

Rq and/or < R^. (2.21) 

For large values of the spectator modulus -R^, the appropriate effective field theory 
description is the {d + l)-dimensional theory of gravity. The modulus Rd becomes the 
Qdd component of the string frame metric, (jdd = (2vr/2rf)^. 

All contributions of \md\ in fc^^ and/or fcy^ are substantial, and the behavior of p is 



better understood in terms of its second expression in (2.11). In addition, and Qy^ 
are exponentially small. In particular, for i?^ 3> -Rg and Rq (which are both > ifl 
one has 



(Coe 
nTf^^^'\z) + nyf^^-''\z) + e''^^+^-\^\\nT + ny)^^ , (2.22) 

where the term e'^*'^^''"''''-' = (Rq/R^^)^ is power-like subdominant. In this expression, 
we neglect terms that are exponentially small in (R^^/Rq) or (R'^^/Rg). The appearance 
of the functions f^~^^^ and fy~^^^ confirms that it is more natural to consider the system 
in (i + 1 dimensions. This is the case since, in this limit, the circle of radius Rd is very 
large. In Fig. [ij the exponential growth of — P when ( — )■ +oo is decreasing. This is 
always the case when ny > but for ny < is only true when z i.e. M/T is small 



enough. This can be seen by considering the large \z\ limit of (2.22) 

-P ~ -T'^el^l-''nTe-2^5^^.2 ^ -oo , 



T"^ e\^\-'^ ny e'^' 82+2 ^ sign(-ny)oo, ^'^''^^^ 

+00 z—>+oo 



-oo 



where S^_^2 is defined in (2.20). We will see that the attraction to the RDS'^"'"^ implies 
z to evolve such that the potential for ( ends by being exponentially decreasing. The 
RDS'^^^ will then correspond to a run away behavior Rd(t) — )■ +oo. The RDS*^"*"^ is 
stable only when the subdominant term {Rq/R'^^Y can be neglected. This is always 
true when Rd is macroscopic, since we restrict our study to temperatures above the 
electroweak scale i.e. Rq not very large compared to the internal space and Rg in 
particular. 



^We introduce e — sign(^) which is 1 in phase III and —1 in the T-dual phase V. 
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If Rd is internal and the subdominant term is not negligible, we will find that the 
universe is attracted back to phase II, where Rd becomes static and the evolution 
becomes that of an RDS*^. Finally we note that at early times, where T and M are well 
above the electro- weak scale, it is possible to have Rd internal while keeping {Ro/R^^Y 
negligible for a large amount of time. In this case at early times the evolution is initially 
attracted to an RDS'^+^. However, at late times {Rq/R'^^Y will always become relevant 
and the evolution always ends in an RDS'^. 

• IV: T-dual fiat potential phase 

— and — < — - Rd, 

Ko rig Kd 

and ^ < Ro and Rg. (2.24) 
Rd 

The effective theory description is T-dual to the phase II where the light degrees of 
freedom are the winding modes instead of the Kaluza-Klein momentum modes of the 
phase II. 

• V: T-duai higher-dimensionai phase 

Ro and/or Rg < (2.25) 

Rd 

This phase is the T-dual of phase III. The light degrees of freedom are the winding 
modes. —P has the same form as in case III, after one transforms ( — > The effective 
field theory is naturally described in d + 1 space-time dimensions, with gdd = i^n/Rd)'^ 
in string frame. 

We have to stress here that the above phases arise from a common string setting but 
cannot be described in the context of a single field theory, due to the lack of the string winding 
modes. In field theory, only the phases II and III (or their T-dual IV and V) can be described 
by a single field theory. In phase I, the string winding modes play a particularly important 
role for the stabilization of the radius Rd at the extended symmetry point, as shown in Sect. 
|3.1 In contrast to field theory, string theory naturally interpolates between these various 



phases. 
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2.3 The type II field theory phases 



The perturbative type II structure of —P can be derived from the heterotic one by taking 
fiT = ny = 0, so that phase I is now equivalent to phases II and III. The local minimum 
of —P at C = is not present anymore and there is a single plateau I U II U IV (see Fig. 
[T]). In phase III (or V), when i?^ ^ Rg and Rq (which are both ^ 1), the function p in 



type II is identical to the heterotic one given in Eq. (2.22). Thus, in type II, the field ( 
admits flat potential phases II and IV in (i-dimensions, and the higher-dimensional phases 
III and V in c? + 1 dimensions. Again the higher-dimensional phases III and V are stable 



only to the extent that we can neglect the {Rq/R'^^Y contribution in Eq. (2.22), which is 
always valid for macroscopic values of i?^. However, when i?^ is internal, the final evolution 
is always attracted to the RDS"^ of phase II or IV. In addition, during the times for which 
{Ro/RdY niay be ignored, the universe admits an earlier evolution well approximated by an 
RDS'^^^. Since there is no enhancement of f/(l) — )• SU{2), the "S't/(2) Higgs phase" does 
not exist perturbatively for type II theories. However, by heterotic-type II duality, we expect 
non-perturbative effects which may enhance the U{1) — t- SU{2) and imply an "5'f/(2) phase" 
I. The non-perturbative effects can correspond to the addition of branes whose separation 
is governed by the spectator modulus. This will be discussed in more details in Sect. |4j 
Alternatively, branes wrapped on a vanishing cycle whose size is fixed by the spectator 
modulus provide another dual type II set up. 

3 Radiation-like dominated solutions (RDS) of the 
Universe and stabilization of the spectator moduli 



In section 2.2, five distinct phases for the thermal effective potential —P were identified 
for the heterotic string. Here, we analyze in detail the behavior of the system in the first 
three phases. The behavior of phases IV and V is found from that of phases II and HI by 
T-duality Rd ^ l/Rd and we do not consider them explicitly. We show that the radius Rd 
can be constant, either at the minimum of the potential in phase I or at any value along the 
flat region II. In phase HI, Rd initially increases along with the expansion of the space-time. 
When the quantities {R^/RdY and {Rg/RdY can be neglected, this evolution continues and 
is well described by an RDS*^"^^. For values of Rd which are internal, Rd is always caught by 
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Ro and Rg and after which the evolution is attracted back to phase II. A distinct RDS*^ exists 
in regions I, II and IV, while in region III and V there exists an RDS"^^^ for macroscopic 
values of i?^. 

Next we show that these cosmological evolutions are stable against small perturbations. 
In particular, for phase I, the spectator modulus R^ is stabilized at the self-dual point, while 
for phases III (and V) it becomes part of the space-time metric. For phases II (and IV) the 
spectator modulus is weakly stabilized due to the presence of gravitational friction arising 
from the expansion of the universe. From these results, one expects in general that when the 
dynamics of all spectator moduli are taken into account, the radii that are not dynamically 
decompactified are (weakly) stabilized at scales smaller than the ones characterizing the 
temperature and supersymmetry breaking. 



For the type II string case described in Sect. 2.3, there is no phase I, due to the absence 
of the heterotic U{1) — SU{2) enhancement at the self dual point. The remaining type II 
phases are identical to the heterotic ones. 

3.1 Case I: Higgs phase 

In this case, the radius Rd is naturally interpreted as a scalar Higgs field for an SU{2) gauge 
group coupled to gravity in d dimensions. As mentioned before, our analysis is restricted to 
field configurations which are isotropic and homogeneous. We thus look for extrema of the 



action (2.3) whose metric, temperature and scalars satisfy the ansatz 

ds^ = -dt^ + aitf(^{dx'Y + ---+{dx'^-'Y^, T{t), zit), (j)±it), Cit)- (3.1) 
The Einstein equations involve a thermal energy-momentum tensor whose components are 



the energy density p and pressure P. Using (2.10) and (B.ll), the energy density p takes a 
factorized form 

p = T'^ r{z,r],() with r = (c? — 1) p — p^, (3.2) 

where Pz denotes the partial derivative with respect to z. Given solutions to the scalar 
equations of motion, we may always find corresponding solutions to the Einstein equations. 



We therefore focus first on solving the scalar equations. Their reduction on the ansatz (3.1) 
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is given in Eqs. (B.35), (B.38) and (B.39) and summarized here as 



oo ^ ° 1 

hd < + (r - p)(f)_i ^ , Prj = 



hC + 



d-2 
1 



d-2 



o 

{r -p)C -Pc = 0, 



(3.3) 
(3.4) 

(3.5) 



where h is defined in Eq. (B.33) and ^ is a function which vanishes when all first and second 



derivatives in its arguments vanish. We have reparameterized our fields in terms of the scale 

o 

factor In a so that time-derivatives have been replaced with (In a)-derivatives denoted as /. 
3.1.1 Radiation-like dominated solution 



To start off, we note that the fact the model is invariant under the T-duality Rd — t- l/Rd 
implies p{z, rj, Q is an even function of so that the first derivative of p with respect to ( 



vanishes at C = 0. Thus, C = is a solution to Eq. (3.5). Next, from Eq. (2.16), the source 



P is independent of at C = 0, so that Prj{z, rj, 0) = 0. As a consequence, Eq. (3.4) is solved 
for any constant = (f)±Q, and we find that (f)± remains a modulus. 

It is convenient to introduce the quantities p{z) and r{z), which are related to the pressure 
P and energy density p at C = as 

P{z,ri,0) =T'^p{z,ri,0) =T'^p{z) p{z,T],0) = T'^riz,r],0) = T'^f{z). (3.6) 



Eq. (3.3) implies the complex structure z can be a constant, long root 



of Vz. As in Eq. (B.40), here Vz takes the simple form 

Vziz) ■.= Vziz,ri,0) = 



d-2 



d- 1 



[r — dp) . 



(3.7) 



The shape of Viz) depends drastically on the model dependent parameter ^ G [—1, 1] 

riT + riT 

and can be inferred from the behavior for large positive or large negative z. 



V{z) 



~ — e 

2— > — OO 



{d-l)z 



1 nv + nv \ , , , 
H — :— I X (Ut + Ut) 



~ -e'^\nv + fiv) X ( 1 

2— > + 00 



2'^ -I riT + fiT 

1 
d 



d 



2{d-l)\l d-1 



d-2 



d-2 

rf-l^^+i' 



with defined in Eq. (2.20). Three cases arise: 
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• Case (a): For ^ — ^ < ^ , V" increases monotonically. 

• Case (b): For < ^ ~^ < 0, has a unique minimum 5c, and p(zc) > 0. 

z'* — 1 riT + riT 

• Case (c): For < ^ — ^, decreases monotonically. 

+ ut 

We choose to concentrate on models where z can be stabilized. This corresponds to the Case 



(b)]^so that 

Case(b): --^ < !^^1±^ < 0, (3.9) 



which guarantees the possibility to fix z at the critical value Zc such that (see Eq. (3.7)), 



f (ze) = dp{Q. (3.10) 

This is the state equation for radiation in + 1 dimension^ The scalar equations of motion 
have now been satisfied and the remaining Einstein equations are easily solved. The overall 



dependence on time is determined from the Friedmann equation (B.33) which takes the form 



hd-2){d-l)H' = ^^ where c, = ^^-^^ r{Q e-"'^ {aoM^Y , (3.11) 
and is easily integrated. Here, Sq and Mq are integration constants. Using the remaining 



Eq. (B.28), the full solution we find is, 

T(t) = M(t) e-'' = e-'-a„M„ , Mt) = <t>io . C(«) = 

a[t) 



(3.12) 



where $ is related to M by the definition (2.8). This particular evolution is characterized 



by a temperature T, a spontaneous supersymmetry breaking scale M, and an inverse scale 
factor that are proportional for all times. From the form of the Friedmann equation, Eq. 



(3.11), this solution is to be interpreted as a radiation-like era in d dimensions, with frozen 



^The models in Case (c) admit a so-called "Moduli Dominated Solution" corresponding to a contracting 
Universe (where z{t) — !■ +oo is running away) [4]. The models in Case (a) are analyzed in [18| and admit 
an RDS in d + 1 dimensions, after dynamical decompactification of the internal radius Rg involved in the 
supersymmetry breaking (i.e. z ^ — 1). 

*As explained in [2j[3], once taking into account the kinetic energy density of the scalar $, one recovers 
the state equation for radiation in d dimensions, as expected for an RDS'^. 
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internal radius Rd = 1 and a modulus (j)±. In the next sub-section, our aim is to analyze the 
stability of this heterotic solution and show that it is an attractor of the dynamics for an 
open set of generic initial boundary conditions. 



3.1.2 Attraction to the radiation-like era and spectator modulus stabilization 



To analyze the stability of the radiation-like era (3.12), we consider small fluctuations around 
it, 

z = Zc + e(,), 0± = 0±o + £^(</,x) , C = + 6(0, (3-13) 

where k(<^x)l k{C)l ^ "^^^ equations of motion for the scalars given in Eqs. 

(B.35), (B.38) and (B.39) become at first order, 

rz - dpz 



°4) + ^(rf-2)l(,)+C'£(,) = where ^=^(^-2)2(^+1) 



OO / 7 \ ° 



0. 



(3.14) 
(3.15) 

2(rf-l) pi~z.) ■ ^^-^^^ 
It is important to note that even if p,^ vanishes and p is independent of rj when C = 0, this 
is not the case for p^^, and indeed we find 



°^°(C) + \id-2) + E{r]) £(0 = where E{r]) 



{d-2f Pa{zc,V,0) 



(3.17) 



Due to the friction term — the solutions of Eq. (3.14) for arbitrary IBC converge to 



as t +00 (with eventually damped oscillations), if and only if C* > 0. This is precisely 



the case when the condition (3.9) is satisfied. Similarly, all solutions to Eq. (3.15) also 



converge to 0. Finally, the generic solution of Eq. (|3.16|) is 

'd-1 



-■iO 



a 



C_L Jd-1 



d — 2 \ai 



d-2 - 
d-1 



.d-1 

+ C_ J d-i , , 

\ d-2 



(3.18) 



where Ja are Bessel functions of the first kind, while c+, c_ (and ai) are constants determined 
by the IBC. As t — )■ +oo or equivalently a — )■ +oo, the above solution converges to with 
damped oscillations, due to the behavior of the Bessel functions. We thus conclude that the 



spectator modulus is stabilized at the self-dual point and that the radiation-like era (3.12) 
is a local attractor. 
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3.2 Case II: Flat potential phase 



We now analyze the phase II. In this case, C is a flat direction of the thermal effective 



potential —P which takes the simple form given in Eq. (2.18). There is no enhancement of 



the gauge group to SU{2) and the spectator modulus Rd is simply a flat direction. Although 
there is no potential for ( in this phase, fluctuations in ( are still suppressed due to the 
gravitational friction caused by the expansion of the universe. In addition, we are going to 
see that depending on the IBC, even if the evolution starts in phase II, it may exit from it 
and enter either phase I or phase III. This result is non-trivial due to the fact that as the 
universe expands, so does the size of the plateau of the flat potential phase. 

As in Sect. |3.1 we restrict our analysis to isotropic and homogeneous d- dimensional 



universes. We take the same ansatz as in (3.1), and consequently all of the equations of 



motion of appendix B.l are vahd in the case considered here. On the plateau, p{z,ri,() is 



equal to given in Eq. (2.18 ), where the dependencies on rj and C are exponentially small 



in Rq or Rg and are thus neglected. As a consequence, P({z, r], Co) — Prjiz, t], (q) ~ 0, for any 
Co on the plateau. 

3.2.1 Radiation-like dominated solution 

Up to the replacement ({t) = — )► ((t) = (q, the vanishing of P(^{z,ri,(o) and pri{z,riXo) 



were the only ingredients used at the beginning of Sect. 3.1 to derive the radiation-like 



dominated solution in d dimensions (3.12). Thus, in the present case a similar solution 



exists. It is obtained by defining p and f as p and r, but with ut and ny set to zero 



P P\hT=nv=0 



r = r f,. 



nT=nv=0 • 



(3.19) 



The condition similar to (3.9) for the existence of a solution with stabilized z is now, 

1 



Case (h) : 



< — < 0, 



2^^ - 1 riT 



(3.20) 



in which case there is a unique Zc satisfying 



f(ic) = dp{zc) . 



(3.21) 
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In phase II, the radiation-like era can be written as 



a{t) 



to) 



ao 



where 



•2 Id 



t, 



<P Cr 



2{d-2){d-r 



1/d 



(3.22) 



T(t) = M(t) e" 



a{t) 



it) 



^±0 



m = Co 



where and Mq are integration constants, is defined as in Eq. (3.11) but with "hat" 



quantities and $ is related to M by the definition (2.8). 



3.2.2 Attraction to the radiation-like era with constant 



The study of the stabihty given in Sect. |3.1.2| can also be apphed in the present case. More 
exphcitly, the equations of motion for the perturbations of the scalars become 



oo 

£( 



) + -(d-2) 



= 0, 
0, 



where C 



{d-2)\d+l) 



Tz - dp^ 



dr — r. 



(3.23) 

(3.24) 
(3.25) 



where the analogue of E{vi) in Eq. (3.16) vanishes, since P((;{zc,'r]Xo) — 0. The positivity of 



C is again guaranteed by the condition (3.20) and all perturbations £(2), £{<j)^), £{q converge 
to as t — 7- +00. 

We conclude that if ( is initially in the range where its potential is fiat, some generic IBC 
imply an attraction to a regime where ( and 0_l are constant moduli. Even if the potential of 
the spectator modulus Rd and are fiat, these fields are nevertheless "stabilized in a weak 
sense" since the expansion of the universe dilutes their kinetic energies which then vanish 
for late times. Physically, the friction terms proportional to H freeze them in place at late 
times. After this has occurred, Rd(t) is a constant, while Rg(t) and -Ro(''^) behave as and 
are thus increasing such that the complex structure-like ratio = Ro{t) / Rg{t) is stabilized. 



3.2.3 FaUing off the plateau 

We now show that despite the expansion of the plateau and the friction terms, for sufficiently 
large initial velocities, it is always possible for ( to escape from the plateau and enter regions 
III or I. To show this, we will simply find particular IBC that yield this precise behavior. 
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Suppose C is somewhere along the plateau and choose, as in Sect. |3.2.1[ z = Zc and 



which are trivial solutions to Eqs (3.3) and (3.4). Eq. (B.30) gives $ 



so that the definition (2.4) implies 



o o , 



T] = {t] + z) 



d-2 
d-l 



d-l' 



(3.26) 



which is nothing but the "constant velocity" of the right edge of the plateau. This defines 

o 

an escape velocity and we now wish to know if it is possible for C, to be larger than this speed 
for a long enough time to reach the edge of the plateau. One can see that taking C, as, 



C = C, + [d-2) 



d 



d 



In a . 



(3.27) 



solves the ^-equation of motion ( 3.5| ) We see that ( is rolling at approximately (see the 



previous footnote) "constant velocity" given by C — — and in particular we have 

o Q 

> Tj. As a result, for initial values Cinit, Vi^t and aj^it, there is always a scale factor Of^n 
where C, reaches the right boundary of the plateau. It is defined by 

IT d-2 



min(?7i„it, t]-,^-,, + z^) - Cmt = \ {d - 2) 



d 



d 



(In a^n - In Oi^iJ > . (3.28) 



Another solution with opposite velocity C, is also allowed, so that C, can roll and enter 
phase I. 



3.3 Case III: Higher-dimensional phase 



Here we analyze phase III of Sect. |2.2[ When Rd is comparable to the size of T*^ ^, it 



is natural to interpret the model in + 1 dimensions. The fields with non-trivial vacuum 
expectation value are the Einstein frame metric g'^^^ and dilaton in d + 1 dimensions. 



coupled to the scalar rj defined in (2.5). We use primes here to denote that the fields are 



now normalized in d + 1 dimensions and are not the same as the fields appearing in the two 
previous sections. Note that Q = \nRct, is now one of the degrees of freedom of g'^^^. In the 
notations of appendix B.2, we are now in the case n = 1 and A = 8 — A — (where A = Q 



for toroidal models and A = 4 for orbifold ones), so that the effective space-time dimension 
isV = d+1. 



^One may worry that this solution yields a singular H . To remedy t his, on e can introduce a perturbation 



HO 



to C, such that |caf^u^| <€. 1, where Ofaii is defined in Eq. (3.28) 
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The action 5* in (2.3) can be written equivalently as 

R' 1 



S = j d^^'xy^ 



2 2^5'^')' -2^^'^'^)' + ^')' 



(3.29) 



where we have introduced the canonically normahzed fields 

2 Id-l 



d 



V 



2 



(3.30) 



v/rfp^l) 

The pressure P' defined in d+1 dimensions is the free energy density —J-"' (see Eq. (B.12)) 
It is related to the partition function Z by, 

Z 



2(d+l) 

P' = e d=T^*dii 



Vr 



0,...,d 



(3.31) 



where the partition function is given as before in Eq. (2.7). In c? + 1 dimensions, the 



temperature T' and supersymmetry breaking scale M' (both measured in Einstein frame. 



see Eqs (B.IO) and (B.5)) have different normalizations compared to their counterparts in 



Eqs (|2.8|), but the complex structure z remains the same. 



r 



e 'i-i 
27ri?n 



M' 



e ''-I 

27lRq 



271 



M' _ Ro 
~V~Rm 



(3.32) 



We are interested in extrema of S that can be interpreted as homogeneous but anisotropic 
space-times with rotation group SO{d — 1). We consider the ansatz 

ds'^ = -de + a'{tf{{dx^f + ---+{dx''-^f^+h{tf{dx''f, T'{t), ^'{t), 0l(t), 

(3.33) 

where the metric scale factor b along the direction d is related to C b> 



b := e'^' . 



(3.34) 



The thermal part of the energy-momentum tensor involves an energy density p', the pressure 
P' (associated to the directions 1, . . . ,d — 1) and a pressure P' + b{dP' /db) in the direction 



d (see Eqs (B.ll )-(B.13)). They are more conveniently written as functions of {T' , z^rj, (} 
in terms of which they take a factorized form, 

P' ^ T"^^'p'iz,r^, Id) , p' ^ T"'^'r'iz,r^, |C|) . (3.35) 

^We introduce absolute values |C| = eC for our analysis to also be valid in the T-dual phase V. 
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The functions p' and r' are related to their counterparts in d dimensions p and r (Eqs. ( |2.10 ) 



and (3.2)) as 



p'(z,r/,|C|) = e''-l^l+^p(;^,r7,|CI), r'(z, r^, |C|) = e''-l^l+^ r(z, r^, |C|), r' = dp' - p'^. (3.36) 



The equations of motion reduced on the ansatz (3.33) are derived in appendix B.2 There 



are three independent Einstein equations, coupled to two scalar equations. It is relevant to 
introduce a modulus ^ as 

■■=-,, (3.37) 
a 

which is a "complex structure for the external space". Then, one can replace one of the 
Einstein equations by the equation for the scalar ^. Effectively, we have three scalars, whose 



equations of motion (B.59), (B.61 ) and (B.62) can be solved, before considering the remaining 



Einstein equations. The scalar equations of motion are summarized here as 



o o ^„ oo oo 



6;'(z,0l,C;^,0i,C;^,0i,C) + K = O 

1 . , , , X °, 1 



h 



oo 



+ 



1 



(p;+pki) = 0: 



(3.38) 
(3.39) 

(3.40) 



where h' is defined in Eq. (B.57) and Q' is a function which vanishes when all its first and 



second derivatives in its arguments vanish. 



3.3.1 Radiation-like dominated solution in d + 1 dimensions 



Our aim is to study the dynamics when the characteristic size of the direction d is larger 



than the scale of the internal space. When ^ Rg and Rq, one observes from Eqs. (3.36) 



and (2.22) that if we neglect the subdominant term e'^^^^^ = (Rq/R^^)'^, the pressure P' 
is independent of r] and ( 



p'{z,r],\C\)-nTf^^'^^\z)+nvfl,^^''{z):=p'{z) when |C| > and r/ + ^ . (3.41) 

In this regime, we define similarly r'{z) := r' and observe that since p'^ ~ ~ 0, constant 
0i = 0io and ^ = ^0 solve trivially Eqs (3.39) and (3.40). 0iand ^ are thus moduh in this 



limit. Note that neglecting the power-like term {Rq/R'^^Y is justified if Rd is macroscopic. 
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since we restrict our analysis to temperatures above the electroweak scale. In particular, if 
Rd takes values such that {Rd/RoY ~ e^'"'^" and {Rd/RoY > e'^'^^^ it is of the order of terms 
which we have already dropped. 



A constant z = z'^ is allowed by Eq. (3.38) if = 0. In the regime we focus on with 



Rd ^ Rg,Ro, the definition (B.60) simplifies to 



v:{z,v,\c\) 



d 



d 



f'-(rf+i)p' ■.= v:{z) 



(3.42) 



which is identical to (3.7) in + 1 dimensions. V'{z) admits a critical point in what we call 
Case (h' ) (by analogy with (3.20)), defined by the condition 

1 



Case (b'J : 



2'^+! - 1 ^ ^ 



(3.43) 



When this is satisfied, z'^ is the unique solution to 

f'iz',) = id+l)p'iz'J, (3.44) 
which is the state equation of radiation in d + 2 dimension^ 

The scalar equations of motion Eqs. (3.38), (3.39) and ( |3.40 ) have now been satisfied. 
The remaining Einstein equations are easily integrated. The overall time dependence is 



determined by the Friedmann equation (B.57), which becomes 

d^ 



-(d-l)dH'^ = 
2 ^ ^ a"^+i 



where 



{d+l){d-l) 



r (4) e 



{d'^r^ (3.45) 



The full solution we obtain is. 



a'{t) 



t \ ''+1 



T\t) = M'{t) e' 



a'it) 



where 



2 

p d+1 



(rf+_l)^^^ 

2{d-l)d 



(3.46) 



hit) 



±0 ' 



where $' is related to M' by definition (3.32), and Oq, Mq, and 0'j_q are arbitrary integration 
constants. For this particular evolution, it is useful to rescale the space coordinate x'^ to 



bring the metric (3.33) in an isotropic form, 

: -de + a!{tY [{dx^Y + ■ ■ ■ + [dx''~^f + [dx"")^ ) . (3.47) 



Id 



ds 



^One has to take into account the classical kinetic energy part of the stress tensor to recover the equation 
of state for radiation in d + 1 dimensions. 
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This shows that there is an enhancement of the local rotation group, SO{d — 1) — t- SO{d). 



We learn from the effective Friedmann equation (3.45) that this evolution of the universe 
can be interpreted as a radiation-like era in d + 1 dimensions. The temperature T' and 
supersymmetry breaking scale M' are inversely proportional to the isotropic scale factor a', 
while (f)'j_ remains a modulus. Next, we study the stability and attraction properties of this 
solution. 



3.3.2 Attraction to the radiation- like era in d + 1 dimensions via decompactification of 
S\Rd) 



To study the stability of the solution (3.46) (which is valid when the subdominant term 



(Rq/R^^) is negligible in Eq. (2.22)), we analyze the behavior of the small fluctuations 
around it, 

z = z'^ + e(^) , 0'^ = 0'^o + , ^ = ^0 + ^(0 , (3-48) 
where )| and \e(^)\ are ^ 1. The scalar equations of motion, Eqs (B.59), (B.61) 



and (B.62), in this regime become 



oo 1 , , O A,, 



oo 



(<^'x) 



^(0 + 2(^-^)^(5) 



where 



C 



^-(d-l?(d + 2) '-'^-(^+1)^ 



(3.49) 
(3.50) 
(3.51) 

(3.52) 



We find that C' always satisfies C" > in Case (b') defined in (3.43). This implies that for 



arbitrary IBC, the solutions of Eq. (3.49 ) converge to (with eventually damped oscillations) 



as t — 7- +00. In addition, all solutions to Eqs (3.50) and (3.51) also converge to 0, when 
t +00. 



Thus, the radiation-like era in d+1 dimensions (3.46) is stable under small fluctuations 



when the subdominant term (Rq/R^^) the Eq. (2.22) is neglected. There is an open set 



of IBC, in particular when R^ is initially macroscopic, such that the solutions are attracted 
by this evolution, which is characterized by an enhanced local rotation group. The Kahler 
modulus Rd is better understood in terms of the "external space complex structure" ratio 
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= h/a = Rd/R, where R is the radius of one of the d — 1 space-hke dimensions. This 
is similar to what is happening to the Kahler modulus Rg that we consider through the 
complex structure-like ratio = M/T = Rq/Rq. However, while z is dynamically stabilized 
to the value z'^, ^ is a modulus ^o- This last remark is due to the fact that we consider local 
equations of motion only. An additional choice of global boundary conditions on the relative 
sizes of the large external dimensions would specify ^o- Only astrophysical observations 
may be sensitive to moduli such as .^o, but not "local" experiments encountered in particle 
physics. 



3.3.3 Residual force and attraction back to phase II 



Here, we want to study the effect of the subdominant term in Eq. (2.22). We take it small 



compared to 1 but not negligible and in particular, we want to analyze how it perturbs the 
results of the previous subsection, namely the attraction to the RDS*^"^^. We first introduce 
y as 



,^+'?-ICI 



Rq 



(3.53) 



It will be convenient to express the conservation of the energy- momentum tensor ( |B.52 ) in 
the following form. 



r' + p') X est. 



Let us consider a generic perturbation around the RDS 



d+l 



Z = K + ^{z) 



^± = 0±o + ' ^ = ^0 + ^(0 , y = yo + £iy), 



(3.54) 



(3.55) 



where \e 



e/y)\ and the constant e'^^" are ^ 1. At order one, the equations 



for £(2) and £{^'^) are identical to (3.49) and (3.50), when e'^^" was neglected. However, the 
equation for e(^^) becomes 



+ ^ (d - 1) ^£(5) + ci e'^^" j = where ci = {ut + ny 



id - l)Sf 

^ ^ >0, (3.56) 



while Eq. (3.54) gives 



e{y) = -C2 £(5) + C3 e(z) where C2 = 1 - 



d 



C3 



4p'(z^) 
1 f^+pl 



d r' + p' 



(3.57) 



This implies e^z) = ^((/>Y) = O5 ^(C) = ~ci e*^^" and = ci C2 e'*^" > 0. We conclude that 
= i?o/-Rd and e^~^ = i^g/i?^ are slowly increasing i.e. the distance between i?^ and the 
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plateau decreases. The interpretation of the evolution in terms of RDS'^"''^ remains valid until 
e'^y and e'^'-^"^-' cease to be <^ 1. As said in the previous subsection, this cannot happen if 
is initially macroscopic and we are interested in temperature and supersymmetry breaking 
scales larger than the electroweak scalej^ However, if initially i?^ is larger than Rq and Rg 
but still microscopic, the above analysis suggests that it will be "caught by the plateau" and 
the cosmology is attracted back to phase II i.e. the RDS'^. 



4 Non-perturbative cosmologies in type II 



As discussed in section 2.3, the thermal effective potential in type II theories does not give 
rise to a Higgs phase I as in the heterotic string. The reason for this difference is that at the 
self-dual point, there is no enhancement U{1) — ?■ SU{2) and thus no growth of the number 
of massless degrees of freedom. However, we expect by heterotic-type II duality that such a 
phase should be possible in type II at the non-perturbative level. A natural candidate setup 
to produce this effect in type II is to introduce a pair of D-branes, whose separation is related 
to the spectator modulus R^. The stabilization of R^ at the self-dual point in the heterotic 
case suggests in the dual type II picture that the effect of the thermal effective potential is 
to fix the D-branes on top of each other, thus producing an U{1) — SU{2) enhancement. 
This attractive force between the D-branes will only be local, in the sense that if we separate 



the D-branes from each other so that R^ is in the range of phase II of Sect. 2.2, the thermal 
effective potential should allow stable finite distances between the D-branes. If we further 
increase the separation, we expect to reach a point where the thermal effective potential 
induces a repulsive force that pushes the D-branes away from each other and we are entering 
phase HI. However, -Ro(^) and Rgit) are increasing faster than Rdit). When they catch it, the 
force between the D-branes vanishes and the distance between them becomes static and we 
are back in phase II. On the contrary, if the distance between the D-branes is macroscopic, 
the effective potential induces a repulsive force that pushes the D-branes away from each 
other and we are in the higher-dimensional phase HI. 

Another set up dual to the heterotic gauge group enhancement can be considered in 



^In non-realistic models (such as the ones considered in this work) where no new physics arises at the 
electroweak scale, considering macroscopic but finite Rq and -Rg catching Rd would correspond to a phase 
of the universe in c? + 2 dimensions with the temperature effectively zero. 
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terms of singularities in the internal space. For instance, a type IIA D2-brane wrapped 
on a vanishing CP^ cycle of radius dual to Rd can give rise to an SU{2) gauge theory and 
admits a mirror description in type IIB [26]. The equivalence between the brane- world and 



geometrical singularity pictures can be analyzed along the lines of Ref. 27 



5 Conclusion and discussion 

In this work, we have considered string theory models in flat space, where geometrical 
fluxes induce a spontaneous breaking of supersymmetry and finite temperature. We have 
computed the 1-loop free energy density, which is nothing but the effective potential at 
finite temperature and first order in perturbation theory. It depends on the temperature T, 
the supersymmetry breaking scale M and "spectator moduli" that characterize the internal 
space but are not involved in the breaking of supersymmetry. Our aim was to analyze the 
dynamics of these spectator moduli in the presence of both temperature and supersymmetry 
breaking. 

We have analyzed in many details heterotic and type IIB models where the dynamics of 
only one of the spectator radii, Rd, is taken into account. More precisely, we have considered 
Euclidean backgrounds which are of the form S^{Ro) x T'^"^ x S^{Rd) x Ai x S^{Rq), where 
S^{Rq) and S^{Rg) both contain fluxes. The flux along the Euclidean time cycle S^{Ro) in- 
troduces temperature T oc I/-R0 and the flux along S^{Rg) implies the spontaneous breaking 
of supersymmetry at a scale M oc I/Rq. The torus T'^^^ is very large, while the internal 
manifold Ai is either T^^'^ or T^"'^ x — , with fixed radii close to the string length. 

In heterotic models, we found five distinct phases of the thermal effective potential (see 
Fig. [1]). In phase I, the potential plays a role in confining the spectator modulus by giving 
it an effective mass in addition to the gravitational friction effects. Rd plays the role of a 
Higgs field stabilized at the enhanced gauge symmetry point Rd = 1, where U{1) — )■ SU{2). 
In phase II (or IV), Rd converges to an arbitrary constant. This is simply due to the 
gravitational friction arising from the expansion of the universe. Thus, while the modulus 
may take any value, its excitations always die off as the universe expands. In phase III (or 
V), if Rd (or 1/Rd) is macroscopic, meaning that we may always neglect the sub-dominant 



term in the effective potential (2.22), then Rd increases proportionally to the expansion of 
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the universe. In this case, the dynamics of the modulus is better understood in terms of 
a complex structure characterizing the anisotropy of a {d + 1) -dimensional universe. As 
in phases II and IV, the excitations of this complex structure die off due to gravitational 
friction. If instead Rd is internal, it always enters phase II (or IV). 

The analysis of the type JIB case is qualitatively the same, up to an important difference. 
The heterotic Higgs phase docs not exist, since there is no gauge symmetry enhancement 
at -Rd = 1 in type II supcrstrings, at least in a perturbative approach. However, we expect 
by heterotic-type II duality that such a gauge theory enhancement should occur once taking 
into account non-perturbative effects in type II superstrings. In particular, the modulus 
governing the distance between D-branes or the size of some cycle on which a brane is 
wrapped could play the dual role of the heterotic radius Rd- 

The heterotic picture of phase I naturally generalizes to the case where all spectator 
moduli are allowed to vary. Although we did not consider this case explicitly, one may 
analyze models with non-diagonal tori. We expect the existence of local minima of the 
thermal effective potential at each enhanced symmetry point, as a consequence of the increase 
in the number of light states. The lowest such point should be given by the most symmetric 
point. 

For the backgrounds S^{Ro) xT'^^^ x S^{Rd) xMx S^{Rg), the thermal effective potential 
has a universal form, up to model-dependent integer parameters {nT,nv)- tit is the number 
of massless boson/fermions pairs in the original supersymmetric model i.e. before the fluxes 
are switched on. ny depends on the precise prescription chosen to break supersymmetry. In 
the heterotic phase I, {ut, ny) has to be replaced by {tit + nr, ny + fiy) to account for the 
additional massless states arising at the enhanced symmetry point. Depending on the IBC, 
we found that the dynamics of the Universe can be attracted to either of the five Radiation- 
like Dominated Solutions associated to the five phases, provided the following conditions are 



fulfilled: 



phase I 



1 



ny + fly 



<0, 



< 



riT + fiT 



1 




(5.1) 



phases II U IV 



2<^ - 1 



phases III U V 



1 




2d+i _ 1 
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When some (or all) of these conditions are not satisfied, different histories of the Universe 
are possible. For instance, suppose a model satisfies the second of the above conditions, 
but not the third, with macroscopic and thus initially in phase III. Depending on the 
remaining initial boundary data, we conjecture at least three different late time behaviors 
to arise, which again correspond to Radiation-like Dominated Solutions: 

• A dynamical compactification of the spectator radius Rd to enter phases II, I, or IV. 
The attractor is an RDS'^, where M oc 1/Rg. 

• A dynamical decompactification of the radius Rg that participates in the spontaneous 
breaking of supersymmetry. This mechanism was conjectured in [2] and is shown in an 



explicit example in 18 . The attractor solution is an RDS'^"'"^, where supersymmetry 



is spontaneously broken by thermal effects only (no M). 

• A (non-perturbative) connection to a cousin model with flux in some of the previously 
spectator directions. For instance, for d + 1 = A and flux in two internal directions, say 

8 and 9, the constraint < — < is replaced by the less restrictive one —0.215 < 

15 n-r 

— < (see |3j|). The attractor solution is an RDS'^"'"^, where M oc 1/y/R^Rg. The 
interesting point here is that the solution is stabilized by the spontaneous generation 
of topological flux. 

In this work, we restricted our discussion of orbifold models to cases where — did not 

Z2 



contain flux and had its radii flxed close to the string scale. In the companion paper 18 
we flll this gap and extend the analysis to cases where an orbifold action is non-trivial on 
dynamical circles and whose radii are either participating in the breaking of supersymmetry 
or are spectators. 

We would like to conclude by giving a summary of the results which have been made 
in studying the intermediate cosmological regime within the framework of string theory, as 
developed in l2|-[5] as well as the current paper: 



The first result is the discovery of the Radiation-like Cosmological Solutions {RDS'^). 
These RDS'^ solutions are not the "usual" radiation solutions defined by pthcrmai = 
{d— l)Pthermai, but iustcad Satisfy ptotai = {d — l)Ptotai only after the contribution from 
the coherent evolution of the supersymmetry breaking modulus is included. These 
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solutions, are also consistent in that the evolution of the space-time curvature scales, 
H"^ = {a/a)"^ and if, the dilaton and the evolving radii scales, (^du)^, (pdii,{Ri/RiY, 
Rj/Rj, is such that they decrease as time progresses. Thus if one starts with small 
curvature and small coupling, at later times one is guaranteed to remain within the 
regime of small curvature and small coupling. 

• The second result is the discovery of the "attractor mechanism" which is valid in the 
intermediate cosmological regime (after the Hagedorn transition or alternatively the 
inflation era but before electro- weak symmetry breaking); within this era the RDS'^ 
cosmological solutions are not only stable under small fluctuations but also are the 
only solutions ( "attractors" ) at late cosmological times. Furthermore, thanks to this 
attractor mechanism, most of the Hagedorn exit ambiguities are washed out in later 
cosmological times. 

• The third result, is the fact that it is possible to derive at the string perturbative 
level, (however exact in a'), the full string free energy J^{T, M; fij) as a functional of 
all moduli (including the SUSY breaking moduli T, M, the string coupling constant 
modulus 0dii as well as the "spectator moduli" ///), at least for a certain class of string 
vacua where the spontaneous breaking of supersymmetry is induced by geometrical 
fluxes. 

There are two main directions in going further with this approach. First, one may carry 

out our analysis with four-dimensional heterotic models, whose internal space in the 

presence of fluxes breaks spontaneously A/4 = 1 supersymmetry. Depending on the details 
of the internal space and spontaneous supersymmetry breaking configuration, it is possi- 
ble for an additional scale Q to appear at very late cosmological times. Q is the "infrared 
renormalisation group invariant transmutation scale" induced at the quantum level by the 



radiative corrections of the soft supersymmetry breaking terms at low energies [19,20 . When 
T{t) < Q, the electroweak phase transition takes place, SU{2) x f/(l) — > f/(l)cm- This starts 
to be the case at a time tw and, for t > tw, the supersymmetry breaking scale M is stabilized 
at a value close to Q. In earlier cosmological times where M{t),T{t) ^ Q, the transmuta- 
tion scale is irrelevant and does not modify our analysis. An important consequence of this 
scenario, is a dynamical explanation of the supersymmetry breaking scale. Indeed, extrapo- 
lating the RDS'^^^ up to the low energy regime where T = 0{1 TeV) one finds, (thanks to 
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the attractor mechanism), that the natural value of the supersymmetry breaking scale M{t) 
is naturally small and around the electroweak phase transition, independently of its initial 
value at early cosmological times. 

The above statement is absolutely correct, if one assumes that there are no other mass 
scales created in the infrared, like for instance the various "infrared renormalization group 
invariant scales Q" associated with: (i) Ac of hidden gauge group (s), or (ii) the transmutation 
scale (s) Qh, (induced in the infrared), by the renormalized structure of the "supersymmetry 
breaking terms". The existence of non-trivial dynamical scale(s) Q modify the very late 
cosmological evolution, namely after the electroweak phase transition around T = Tw ~ 
Qh ~ 0{lTeV). Due to this, the intermediate cosmological regime which we study in this 
paper is defined by : Te ^ T(t), M(t) ^ T\y- 

Obviously, the physics for t ^ t^r is of main importance in (astro)particle physics and 
late time cosmology. Unfortunately, the infrared phase at t ^ tw depends strongly on the 
specific choice of four-dimensional = 1 superstring vacuum and the way of spontaneously 
breaking supersymmetry. A lot of work is necessary to determine the initial superstring vacua 
that leads in late-times to the precise structure of our Universe. On the other hand, we would 
like to stress here that the qualitative infrared behavior of the effective stringy "no-scale" 



field theory 20 28 strongly suggests that we are in an interesting string evolutionary scenario 
after the "Hagedorn-transition exit", t > Ie, connecting cosmology to particle physics. 

Secondly, in the regime close to or at the Hagendorn transition, the conventional notions 
of General Relativity such as geometry and topology are well defined only in the low energy 



and/or small curvature approximations of a string theory setup 29 . In the very early times of 
the Universe, t < Ie, purely stringy phenomena at very small distances and strong curvature 
scales imply that the physics could be quite different from what one might expect from a 



"naive" field theory point of view 29 . In this early epoch, classical gravity is no longer 
valid and has to be replaced by a more fundamental singularity-free theory such as (super- 
)string theory |16]. Thus, the main obstruction in such a stringy cosmological framework is 
the Hagedorn temperature limitation T < Th- Actually, this is not a pathology but rather 
the signal that a phase transition from a previous vacuum is taking place. The Hagedorn- 
like singularities have to be resolved either by a stringy phase transition |6|[8}[T2|[T3] or by 
choosing Hagedorn-free string vacua in the early stage of the universe [TslfTG]. 
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In this work, we have bypassed the Hagedorn transition ambiguities by considering arbi- 
trary initial boundary conditions (IBC) at tE, the "Hagedorn transition exit time". Thanks 
to the attraction to "Radiation-hke Universes" in late times, most of the ambiguities are 
washed out. It is however of fundamental interest to investigate further the early non- 
geometric era of our universe and to show that the induced IBC at tE solve naturally the 
"flatness" and "entropy" problem in late cosmological times. This stringy scenario would 
be an alternative (or at least complementary) to the conventional inflationary scenarios 
proposed in field theory. 
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Appendix 



A Partition function 
A.l Heterotic string on tori 

We first focus on tlie lieterotic string compactified on a Euclidean toroidal space, while 
heterotic and type IIB orbifolds will be considered in appendix A.2[ Our starting point 



is the heterotic string in a background x x T" = 9 - D), where S}^{Ro) is the 
Euclidean time compactified on a circle of radius Rq- For simplicity, we choose the tori 
and T^-^ to be products of circles Ilp^^S^{Rp) and Ilf^^_^_lS^ (Ri) . The partition function 
vanishes, due to supersymmetry: 

a,b ' ' p=l i=D+l 

(A.l) 

To implement finite temperature, we deform the i)(_Ro) lattice by coupling the space-time 
fermion number Qp = a to the momentum along S]^{Rq). We also introduce a spontaneous 
supersymmetry breaking by coupling R-symmetry charges a + Qi {i = D + 1, . . . ,9) to the 
momenta along the T" directions. Spin statistics and modular invariance then determines 
the precise replacement of the lattices as 

* m,n 

(A.2) 

/{, is the set of labels associated to directions with fluxes that break spontaneously super- 
symmetry. In practice, the Qj's ii G /;,) are linear combinations of charges of the x E'^ 
lattice, for which Si is determined to be or 1. In our notations, Qo = -^o = and eo = 1- A 



convenient rewriting of the phases in Eq. (A.2) is done by defining mj = 2ki + gi, Ui = 2li + h. 



and summing over Qi, hi G {0, 1} and fcj, Zj over all integers. We may evaluate the sum over 
the spin structures a and b by redefining a = a + J2ieit ^* b = b + J2ieii, 9i- "^^^ phases 
from ( |A.1[ ) and ( |A.2[ ) combine to give 

a + 6+a6+^ (jii{a + Qi) + hi{b + Li) + eiQih^ = d + b+ab+y^Ji{l + Qi) + P{g, h, Q, L, e) 

(A.3) 
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where P{g,h,Q, L,e) consists of terms which vanish when hi = (Vi G /;,). We may now 
make use of the Jacobi identity 

1 ^(_l)a+S+aS^4[a-.g^j ^ -^^Kli^] (A.4) 



to obtain the result 



2 

d'b 



dTidT2 



nr(i,i)S:](^.)xnr(i.i)Wx^, (A.5) 

ieh p=i ' ' 



where we have introduced the shifted lattices 



1(1,1) 5;] (i?,) = ^^Y^e-^i\^''+9M+hM\ (A.6) 



Important simplifications can be made, using the fact that -Rj ^ 1 (i G // 



b ■ 



When 2li + /ij 7^ in Eq. (A. 6), the integrand in the partition function ( |A.5 ) contains 



a factor e '^^i (2'i+'i«)t-2 implying an exponentially suppressed contribution after integration. 
We thus restrict to the sectors with Zj = /ij = (z G /;,). 

- The sectors hi = Qi = mod 2 are supersymmetric and therefore do not contribute, 



as can be seen from the presence of a '&'^[\] factor in (A.5). We thus only keep the sectors 
hi = {i e h), Yji9i = 1 mod 2. 



- All of them, have at least one i E h such that Qi = 1, so that the integrand in (A.5) contains 
a factor e ^2 . This implies that we can extend the integral over the fundamental domain, 
to an integral over the full upper half strip. The error introduced this way is exponentially 
suppressecj^ 



^Note that if the model before introducing temperature and internal fluxes was not supersymmetric, the 
sector hi = Qi = Q (Wi S h) would not vanish. Its integral over the fundamental domain would not be 

suppressed by any e ^2 (i g /j,) factor and could not be replaced by the integral over the strip. The result 
would imply a very large (but finite) contribution to the vacuum energy, proportional to the number of 
massless bosons minus the number of massless fermions. 



37 



Altogether, the partition function reduces to 



dri 



+00 



9t- 2 

T,i 9i = lmod2 



(A.7) 



D 

p=i 



In this expression, the low lying contributions from the oscillators and right moving 
lattice r(o,i6) are 



{k.i 



where q = e'^^^'^ and D^icj^Q) is the sum over massless degrees of freedom with each mode 
weighted by the factor (— l)Efc9fcQfc_ Defining = {1, . . . , D} the set of "spectator" directions 
i. e. that do not break supersymmetry, the lattice of can also be expanded as 

D 

J]r(i,i)(/2p) = e-""^^''i^^^'+^"^''''^')e-2*-^i^.™-"^ (A.9) 

P=l mq,nq{qels) 



The change of variable T2 = xtt RfC^^i + 5'*)) the r2-integration shows that the 



massive contributions in Eqs (A.8) and (A.9) are exponentially suppressed, compared to the 



massless ones. We thus concentrate our attention on the hght degrees of freedom. Because 
we have been able to replace the fundamental domain with the half strip, the integration 
over Ti now simply enforces the level matching condition. The constant term in the r.h.s. of 



Eq. (A.8) combined with Eq. (A.9) implies that it is enough to keep for each p G Is 





rir, 



nip = 



if i?„ ^ max Ri 



1 

if Rp < 

max Hi 



(A.IO) 



nip = rip = else. 



Similarly, the q ^ term in Eq. (A.8) combined with Eq. (A.9) implies that it is enough to 
keep for each q E Is, 



m„ = nr, 



±1 if Rq ^ 1 with mp, rip given in Eq. (A.IO), Vp G J^,]? 7^ g. (A. 11 



Physically, the above two winding modes are responsible for the gauge symmetry enhance- 
ment U{1) —7- SU(2) at Rq = 1. Away from the self dual point, they are massive: Their 
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contribution to Z becomes exponentially negligible, while the SU{2) is Higgsed. The en- 
hancement of the symmetry will play an important role in stabilizing Rg around one. The 
following reduced expression for the partition function when -Rp > 1 (Vp G Is) is then 
obtained, 

n I E I r ^ f-f E + 

DoilQ) E -p('-vrr,5^(^)') 

nirirels) \ pels ^ P/ / 

+2 5:exp(-.r.(i, + fl^2)) ^ exp f-.r. ^ (|^)^ 

(A.12) 

The companion expressions when any i?p < 1 (p G /g) is obtained by replacing i?p — — . 

Rp 

The above result can be rewritten in various forms. We introduce an integer parameter 
A and split the torus as T^^^ x T^. Correspondingly, the "spectator" indices Is are 
divided in two sets, 

= {1, A} , = {/} - A + 1, . . . , D}, (A.13) 

and we perform a Poisson resummation on the T^~^ zero mode lattice. The reason for that 
is that when some Rp {p G Is) is "small" i.e. Rp <^ infjg/j, Ri , the Hamiltonian formulation 



of Eq. (A.12) is relevant, while when some Rp {p G Is) is "large" i.e. Rp ^ maxjg/j^ Ri, the 
Lagrangian formulation is more convenient. The alternative forms for arbitrary A are: 

^=(n^^) n^^' E i^l'!.+3 exp ( -- >^ RK2k + ~gi)' 

9j = lmod2 




E I 7 E ^^p'^p'f ~ E 

mr{r £ l^,r q) 
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TT 



mr(r G if) 



(A.14) 



Using the integral form of the modified Bessel function of the second kind Ka{z) 



(g-l) 

j^^ ^ exp (-^^) exp(-7rr2G) = ' 2ir„_i(27rv^), (A.15) 



one obtains 



^ = n n 



E 



mr(r &ll,T ^ q) 

+2E E 

g'e/| m,,/(r' e /^jr' ^ g') 
m,.(r e /|) 



Ej 9j = Imod 2 



_D-A + n + l 



D-A+n+1 



(27r/F\G^) 



( -D-A+n 



(A.16) 



where we have defined 



1 „ \ sr^ I m.„ 

Kq' 



E \ R 



(A.17) 



In the second hne of Eq. (A.16), the "primed" sum in the brackets means that = 



(Vr G Ig) is excluded. We remind the reader that Eq. (A.16) is valid when Rp > 1 
(Vp G U II). The expressions with some Rg's such that Rg < 1 is obtai ned by T-duality 



z.e. by exchanging them with their inverses, Rg — )■ — , in Eqs (A.16) and (A.17). 

Rn 
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It will be convenient to have the rules for decompactifying directions of ^ as well as 
the rules for freezing one of the radii. 

• For decompactifying a radius Rg' {s' G /'), one simply keeps only the terms with 



rhs' = in (A. 16). In the last line, one also discards the term with q' = s'. The net 



result is a remaining overall factor of Rg' in the first line. 

In order to freeze a radius Rg (s G /|) at the self dual point Rg = 1, one keeps only 



the terms with rUg = in (A. 16). In addition, one discards the term with q = s in the 



third line and shifts Dq ^ Dq + 2. 

In order to freeze a radius Rg {s G /|) at an arbitrary value such that l/Rg and 



Rg <^ inf jg/j^ one keeps only the terms with nig = in (A. 16). In addition, one 
discards the term with g = s in the third line. 

A. 2 Heterotic and type IIB orbifolds 

We next consider the heterotic string on S]^ x x — x (n = 5 — D). The partition 
function vanishes, due to the 8 conserved space-time supercharges. 



2^2^' ' Wl-lb+Gi-lb-Gl-y'^,'^)lUl 8^20 

^ H,G a,b ' ' 

D 9 



(A.18) 



xr(i,i)(i?o) X J]r(i,i)(i?p) X II r(i,i)(i?,), 

p=l i=D+5 



where H, G are equal to or 1. -Z'(4^4)[^] is the block that accounts for the — part of the 
background and the right moving r(o,i6)[G] lattice is consistently Z2-twisted to guarantee 
modular invariance, thereby breaking part of the initial gauge group. Temperature and 
supersymmetry breaking are again introduced by modifying the lattice sums along S]^ and 



T" as in (A.2). However, the left moving charges Qi may now also involve the "orbifold twist 



number" H. 

The analysis of the toroidal case can be applied the same way. Defining = {0,1^ + 
5, . . . , 9}, any sector with some /ij 7^ (i G h) is exponentially suppressed and the sectors 
with Ylieii, hi = and X]je4 9i ~ ^ vanish due to supersymmetry. We may again replace 
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the fundamental domain with the full upper half strip. We treat — as part of the internal 

Z2 

sector, with frozen moduli much smaller than inf As before, the non-exponentially 

suppressed contributions to the partition function arise from the massless modes (and their 
light towers of KK (or winding) states). For explicit examples, see [2]. The net result is that 



the partition function is of the same form as in (A. 16) except that the numbers Do{g, Q) are 
different. 



rp4 

The type JIB partition function on 5']^ 

Z2 



(n 



5 — D) takes the form 



dTidT2 1 1 

2 ' 2 

H,G a,b 



2T2 



-1) 



X 



^(4,4) 



D 



ib-c] 2 



ib-G- 



a,b 



xr(i,i)(i?o)x J]r(i,i)(i?p)x W r(i,i)(i?,). 



i=D+5 



(A.19) 

Temperature and supersymmetry breaking are again introduced by modifying the lattice 
sums along S]^ and T". However, in the type IIB case, we may introduce phases similar to 
Eq. (A. 2) but involving either left moving R-charges a + Qi (z G h), or right moving ones 
d + Qi, or both. In the present paper, we consider cases where both left and right charges are 
non-trivial. Some models in this class where shown to allow critical cosmological evolutions 
corresponding to radiation-like eras [3||4j. 

For Ri ^ 1 (Wi G h), the manipulations used in the heterotic case can be applied 
similarly, up to an important difference. In the sector [q] = [q], the analogue of the heterotic 



contribution given in Eq. (A.8) is in type IIB, 



2'il + Oiq,q) 



(A.20) 



Consequently, there is no massless winding mode arising when some i?g = 1 (g G Is), as 



opposed to the heterotic states given in (A. 11). The final form of the partition function 



is then formally identical to the two first lines of the heterotic one (A. 16), with coefficient 
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Doig,Q,Q), 



\keib 



X ■ 



J2j9j = lmod2 

(ttFi) 2 



D-A+n+l 



-A+n+i (27rA/ -FiGi , 

2 

(A.21) 



where Fi, Gi are defined in Eq. (A. 17). 



B Equations of motion 

In this appendix, we derive the equations of motion in our thermal backgrounds. The 9 
space-like directions are x T" in the toroidal models and x — x T" in the orbifold 
ones. To unify the two cases, we define n = 9 — A — D , where A = in toroidal compactifi- 
cations and A = 4 in orbifold ones. Geometrical fluxes in the T" torus break spontaneously 
supersymmetry. A T'^~^ sub-torus of is taken to be isotropic and very large, to be inter- 
preted as part of the spatial directions of the space-time. In practice, this means that the 
radii of T'^~^ are proportional for all time t and that the associated KK states are continu- 
ous, implying the convergence of the partition function. The remaining 

rpD-d+i sub-torus IS 

allowed to have arbitrary radii. When they are all small, we interpret them as internal and 
the space-time dimension is d. When only A of the T^^'^^^ radii {D — {d — 1) > A > 0) 
are small, the space-time dimension is then V = D + 1 — A. However, this space-time 
with enhanced dimension is anisotropic since only part of its space-like radii are evolving 
proportionally. 

Our starting point is the standard 10-dimensional dilaton-gravity theory, 

S=lj Av^e-2*^"^° [i?io + 49^0cmio5'^0dnio] " j d'^Xy/^oAo- (B.l) 

The hats denote that we are in string frame and the numerical subscripts, here 10, indicate 

the space-time dimension. 0diiio and J^iq are the dilaton and free energy density in 10 

Z 

dimensions. The latter is related to the partition function by J-'io = ~t^, where Vio is the 
10-dimensional volume of the Euclidean background in which we computed Z. We split the 
torus as T^^^ x and dimensionally reduce on and T" (or — xT"" in the orbifold 
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models). The action becomes 



-gve 



-20dilI5 



D 



p=D-A+l P i=D+l 



V, 



(B.2) 



where the dilaton in P-dimensions is 0diii> 



^diiio - I Y^l=D~ A+i ln27ri?^ and = - — ■ 
In some instances, we will suppose for simplicity that some internal radii are frozen. This 

will be the case for some radii of T^, and the radii of — in the orbifold models. Then, 

Z2 



the corresponding kinetic terms would simply disappear from the action (B.2). Performing 



the conformal transformation g-ri = G^p{—j^4>diiv)gv brings us to Einstein frame and the 
action becomes 



S = ^ I (fx^J-gv 



R 



V-2 ^ 



p=D-A+l P i=D+l 



j (Px^-gvTv-, 



(B.3) 



where Tv = exp(^|^</)diiX))-FD. The supersymmetry breaking scale measured in Einstein 
frame, M-p, is given by the inverse volume of T", 

9 



, ^ ^'/'dilP -I— r 

Mv = e~^ [[ 



1 



(27ri?i)Vn 2n 

- 1 



exp 



i=D+A+l ^ ' ^ i=Z)+yl+l 

It is convenient to introduce an explicit field notation, $x), for the supersymmetry breaking 
scale as 



V = where = - — - + -. (B.5) 

The coefficient a is chosen so that has a canonically normalized kinetic term. We can 
introduce other fields, and (i = 1, . . . , n — 1), to describe the remaining degrees of 
freedom among the dilaton and radii of T". The explicit transformation law is given as 



V ^1 / 



V-2 



1 

an 

1 



_ J_ 

an 

1 



__!_ 

an 
1 



V+n-2 y/V+n-2 

"-^ ^ 



-^J n(n— 1) n(n— 1) 





yJV+n-2 
1 

a/"("-1) 



1 

V2 



lni?9 
lni?s 



.(B.6) 



\/2 J 



\ hiRo+A+i J 
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Finally, we denote the A degrees of freedom of as (p = In Rd^a+p (p = 1, • • • , A). In 
terms of these fields, the action takes the canonical formp^j 



n-l , A 



L i=l p=l 

- j (fx^^Fv. (B.7) 

The metric variation gives the standard Einstein equation in terms of the stress energy 
tensor. It is convenient to single out the terms coming from the thermal corrections and so 
we define the thermal energy-momentum tensor as 

dJ'-D 

Tv^u = -gv^^uJ^v + 2-^^—^. (B.8) 

We are interested in space-times which are homogeneous but anisotropic, since the radii of 
rpv-d allowed to vary independently of the radii of the very large T"'"^. In this case, the 
fields are only allowed to vary with time and we take the metric to be of the form 

d-l V-d 

dsl = -dt^ + a{tf ^dxi + ^ hk{tfdxl_^_^i,. (B.9) 

1=1 k=l 

The hk are the scale factors of the T-^"'^ torus in Einstein frame and are related to the radii 

— 2 ■ 

in string frame by hk = ei'-2'^'*"^27ri?fc if -Rfc > 1 (exchange — )■ 1/-Rfc if ^ !)• It is also 
convenient to introduce the temperature related to the radius of the Euclidean time used in 
the computation of Z as 

(not to be confused with the stress energy tensor). Defining the thermal energy density 
p-D = Tdoo and pressure Pv = a^^T-Dn (no sum on / = 1, . . . ,d — 1), the thermal energy- 
momentum tensor can be expressed as 
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r.„o ^ -F. + 2r|§|^(T.^-P.) (B.ll) 

Tvii = -a^J^v = a^Pv (B.12) 
Tvd+kM>^ = -^l^^ + ^^2=bl[Pv + h^y (B.13) 

-'^"One could also introduce A — A more scalars associated to the radii of /'L2 in the orbifold models. 
However, as announced before, we consider from now on these radii to be internal and constant. 



^^The extra temperature factor of in the second term of the first equation in (B.ll) can be seen 
as follows. The variation of the free energy density can be taken when the coordinates are such that the 
metric is just the analytic continuation of the Euclidean background, ds^ = —T^^dt^ + a{t Y ^sj + 



S/c ^fc(*)^'^*^d-i+fe- '^h^ factor of then comes from changing the metric coordinates to have (B.9) 
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Note that the thermal energy density satisfies the thermodynamical identity p-p = 

T-J^—^ — Px). The Ricci tensor can be expressed in terms of if = - and K]^ = — 
oTd a bk 



,V — d). The off-diagonal elements automatically vanish, as well as those of 



the thermal energy- momentum tensor (B.ll). The remaining diagonal Einstein equations 
become 



P - 3 



V-d 



n~l 



i=l 



^ k=l ^ ^ i=l 

k=l ^ ^ 



V-2' 



(B.14) 



k=l 



V-2 



V-2 



1 ^''^ f 



k=\ 



V 



dbk 



(B.15) 



v-d 



V-2' 



V-2 



Pv + b, 



dp 



V 



dbk 



d-1 
V-2 



1 '^ '^ / dP \ 

+ E + {k = l,...,V-d), (B.16) 



where D = + 1 — A. The equations of motion for the scalars reduce to 

v-d 



^v + {d- l)H^v + J2 ^^^^ = 

k=l 
V-d 

+ {d-l)H^^j, + J2Kk^ ±v 

k=l 

V-d 

C^ + {d-l)HCp + Y.KuCp 

k=l 

v-d 

(pi + {d- l)H(pi + ^ Kk^. 



dP-i 



V 



dPv 



■>±v 



dPn 



V 



k=l 



dCp 

dPy 



(p=l,...,A) 
{i = l,...,n- 1). 



(B.17) 
(B.18) 
(B.19) 
(B.20) 



B.l Reduced equations of motion for phases I and II 



Here, we apply the above results to the background of Sects 3A_ and 3^ We have n = 1 
internal direction that breaks supersymmetry and D = 8 — A "spectator" directions {A = 



46 



for the toroidal models and A = A for the orbifold ones). We split = T'^ ^ x (^S^{Rd) x 
T^-'i) , where T'^-^ is very large, the radii of T are small and frozen, and we are interested 
in the regime where Rd and l/Rd are smaller than inii^j^ Ri. We thus choose A = 9—A—d and 
find an effective space-time dimension V = d. Beside the temperature Td, the independent 
fields are the scale factor a and the scalars (p±d and Ci = In Rd, whose expressions follow 



from Eq. (|R6|, with a = ^{d - l)/{d - 2), 



Ci 



V(rf-2)(rf-l)"' 
2 1 



d-2 



Vd^ 
\nRd. 



d- 1 
In Rg 



In Rq 



(B.21] 



To simplify the notations, we will denote the fields as T, a, C and the thermal energy 

density and pressure as p, P. 



The two Einstein equations (B.14) and (B.15) simplify and can be replaced by the Fried- 



mann equation and an equation expressing the conservation of energy, 

^(d-2)(d-l)/f2 = ^($^ + 0i + C')+p, 



p+{d-l)H{p + P) + ^ 



dP ■ dP -dP 



(B.22) 
(B.23) 



where the sources P, p satisfy (see Sect. 2.1) 



P = T''p{z,r^,0, P = T 



dP 

df 



p = T'^r{z,r],C) with r = {d-l)p-p„ (BM) 



where ?7 = Ini^g and is the ratio of temperature to supersymmetry breaking i.e. z 



- ln(27rT) 
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The scalar field equations (B.17)-(B.19) reduce to 



^ + {d-l)H^ 

1 + (rf - l)H<i)^ 
C + {d-l)HC 



9$ [y d-2 

dP T'^ 



Pz 



d-2 



d 



dP 
9C 



Vd^ 



(B.25) 

(B.26) 
(B.27) 



"'^^In this section, it is understood that partial derivatives of p are with respect to z, rj or ^ with the 
remaining variables held constant. 
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, , dP ■ dP ■ dP ■ dP 

Note that Eq. (B.23) can easily be integrated. Writing P = $ — — — — hC -ttt + T — — , 



one derives from Eq! (B.23) and the relation between p and P in Eq. (B.24), 



p + P + (t/-l)i^(p + P) = ^(p + P) 



{aTY-\r{z,r,X)+p{z,nX)) =csi. (B.28) 



It is useful to parameterize our functions in terms of In a and thereby replace time- 
derivatives by (In a)-derivatives, in which case we have 



f = H 



df 
din a 



Hf. 



(B.29) 



for any function /. Critical solutions do not have constant $ but rather constant z and so 



it is relevant to change variables from $ to 2 in Eq. (B.25). Physically, this corresponds 



to the fact that for stable solutions the ratio of the supersymmetry breaking scale to the 
temperature scale must be a constant. In order to proceed, we first express the derivative of 

o o o 

the energy density p in T, z, t] and ( variables as p = T'^H{drT /T + rzZ + rrjrj + r(^() and note 

° / ° ° 

from the definition of z that T/T = ^J{d — l)/{d — 2)$ — z. Using these two expressions, 



Eq. (B.23) can be reexpressed as 



(B.30) 



where 



Vd^ S 

r + p 



8 ' 



and 



(B.31) 



(B.32) 



Eq. (B.30) can now be used to eliminate $ from the equations. We first use it to write the 



Friedmann Eq. (B.22) in the form 



H = hT where /i(z, (; 2, C) = y 



{d - 2){d - 1) - IC 



and 



o o 
,2 I /■2 



(B.33) 



(B.34) 
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Next, noting that $ = if$ + one can express H in terms of H,p,P using Einstein's 

equations and bring (B.25) into the form, 



h 



.00 . 



+ A^QC + il<P±,C)C 



+ 


h 




+ 


h 


( ^ 


+ 


h 


(Be- 



d- 



(B.35) 



d-2 



f (Ac)^)J + ^(r-p)Ac) 



where the matrix C is 



^ Az)^ — \ l %AAiz\A 



A 



{z)r) 



/d-1 



A< 



A(d,,)7. \l j^_^A(z)A(^if,^)ri A{^^^)r] { rj—i 



d-2 

d-l'^('P±. 
d-2 



Ac 



j^A(z)A[(^yr) A(n'„ f — \l ^A( 

















1 ^ Af) Ac)'? 





o ° ° 

and we have introduced V{z, r], (; z, C); whose derivative with respect to z is 
K = • 



d- 1 



d-2 



Pz 



d-2 



d-l^" ' d-2 



^ 'r-p)B- J'^^hBr.B. 



d-2 



(B.36) 



(B.37) 



Similarly, the equations (B.26), (B.27) for 0_l and C, become 



00 ]^ ° 1 



d-2 



Vd^ 



hC + ir-p)C-Pc = 0- 



(B.38) 
(B.39) 



Finally, note that for ( = 0, Vz defined in (B.37) does not depend on velocities anymore. 



B.2 Reduced equations of motion for phase III 



We want to write the fields and equations of motions in the regime of Sect. 3.3, Again, 
n = 1 and D = 8 - A. We split as {T'^-^ x S\Rd)) x T^""^, where T"^"^ is very large, 
the radii of T^~'^ are small and frozen, and we are interested in the regime where Rd (or 
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its inverse) is large. We thus take A = 8 — A — d, which imphes an effective space-time 
dimension V = d + 1. The independent fields are the temperature T^+i, the scale factor 
associated to the torus T'^~^, the scale factor bi of the spatial circle S^{Rd) and ^d+i, (p±d+i- 
The definitions of the scalars is derived from Eq. (B.6), with a = a/ d/{d — 1), 

2 



d+l 



2 1 



d-1 

~~d 



We introduce simpler notations for the fields in d + l dimensions, T', a', b, and for the 

thermal energy density and pressure, p',P'- 



There are three Einstein equations (B.14)-(B.16) which simplify as 

1 



{d - 1){H' + H'^) - (K + K^) 
H' + {d- 1)H'^ + H'K 
k+{d- l)H'K + K"^ 



d- 1 



{d-2)p +dP' + b 



dP' 

~db 



d-1 
1 

d- 1 



p'-P'-b 



dP' 

~db 



p' - F + {d-2)b 



dP' 

~db 



(P.42) 
(B.43) 
(B.44) 



where H' = — and K = -. The dependencies of the thermal sources P' and p' are as follows 
a' b 



(see Sect. 3.3^ , 

P' = r'^+V(z,r^,|C|), p' = T 



, dP' 



P' =^ p = T"^'^^ r {z,'r],\C,\) with r' = dp —p^, 

(B.45) 

where rj = Ini^g, ^ = Ini?;^ is related to the definition of b via |C| = ln6 + — ^==$' + ~^4>'± 
and is the ratio of temperature to supersymmetry breaking i.e. z = 



d ^1 



j<l>'-ln(27rT') 
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The scalar equations of motion given in (B.17) and (B.18) become in the present case, 

dP' 



$' + (rf - 1)H' + K]^' 



{d-l)H' + K) 0'_ 



9$' 

rjnld+l\ 

dP' 



d-1 



d . d — 1 . 

Pz-\l—^Pv + 



1 



d 



^{d-l)d 



501 

rpld+l 

^Td 



(p;+Pici) • 



(B.47) 



^^In this section, it is understood that partial derivatives of p' with respect to z, r\ or \C\ are with the two 
other variables held constant. 
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Again, it is convenient to derive from Eqs. (B.42)-(B.47) the Friedmann equation and the 



conservation of the energy-momentum tensor, 
1 



{d-l){{d- 2)H'^ + 2HK = - $'2 + + P' 



(B.48) 



P' + {(d - 1)H' + k) ip' + + ^' + 01 + & ^ = 0. (B.49) 



(90'^ db 

To obtain a useful form for the remaining independent equations, we first define 

b_ 
a' 



(B.50) 



and then subtract Eq. (B.43) from (B.44) to obtain, 

- / \ ■ BP' 



(B.51) 



• , dP' dP' ■ dP' ■ , dP' 
Eq. (|B.49|) can be integrated by writing P' = $' — — + 0', -— - + h — — h T' —— to obtain, 
^ ' a$' 0(p, do oT' 



p' + P' + {dH' + Oip' + P') = ^ip' + P') 



[a'T'fe^ (r'(^,r/,|C|)+p'(z,r/,|C|)) =cst. 

(B.52) 



As in appendix B.l, we introduce (Ina')-derivatives, where for any function /, 

w df 



f = H' 



din a' 



H'f. 



(B.53) 



Proceeding in an analogous manner as the derivation of (B.30), Eq. (B.49) can also be 
rewritten as 

I' = Al^z + ^'(^^)0i+ A\J + B' , (B.54) 

where 



{d +l)r' - 
S' 



'^'{cf>'^){z^v, ICI) = - 



Vd £' 



(B.55) 



r' + p' + r[.| + p\., , d(r' + p') + r',., + p\^. 



and 



d 



d-1 



8' 



((rf+l)r'+p;) 



8' 



51 



Eq. (B.54) can be used to recast the Friedmann equation (B.48) in the form, 

r' 



where 



IC' 



where i], |C|; 2;, = t 



{d-l)d-lC''' 



(B.57) 



(B.58) 



For the scalar its equation becomes 



h' 

+ 
+ 
+ 



/n 



A' B'] +A' ] + ^'"^"^i^i^ V 



h 



I / 1 



B' + i3' , - 



h! B\ 



ICI 



d-1 



(0 



(B.59) 



where the components of the matrix C are collected below, in Eqs (B.63)-(B.65) and 

o ° ° 

y{z,r], ICI; z,0_L,C) is defined by 



V = - 



rf- 1 



d I d— 1 I 

Pz+\l —^Pv 



d 



^{d-l)d 



PiCI 



+ —{r'-p'-p\^^)B' + h' 



B', 



ICI 



B' + l 



— B' B' 

d " 



(B.60) 



The remaining equations (B.47) for (p'^ and (B.51) for ^ take the form 

1 



00 1 



h' (h\ + 



h 



d-l^''-^'-^'\^\^^'- ^d 
1 



+ ^ (^'-p -pjci) e-pjci = 0. 



(B.61) 
(B.62) 



For completeness, we collect all entries of the matrix C that appears in Eq. (|B.59|): 

1 



C = A!, ^ - A', X 



z - A4>'^)z - Az) 



d-l 



A, X - 



A!< 



d ^^^^ v^(d^^ 



^)ICI 



^ A' 



(B.63) 



iz 
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z'/ _ l3ii!7__Mii 4' W- tAi A' 



= Aon - Ao [^I^Aov - ^(^^•^(OICI I ■ 
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